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We study the antiferromagnetic classical Ising (AFI) model on the sorrel net, a |th site depleted 
and =th bond depleted triangular lattice. Our classical Monte Carlo simulations, verified by exact 
results for small system sizes, show that the AFI model on this corner-shared triangle net (with 
coupling constant Ji) is highly frustrated, with a residual entropy of jL = 0.48185±0. 00008. An- 
ticipating that it may be difficult to achieve perfect bond depletion, we investigate the physics 
originating from turning back on the deleted bonds (J2) to create a lattice of edge-sharing triangles. 
Below a critical temperature which grows linearly with J2 for small J2, we identify the nature of 
the unusual magnetic ordering and present analytic expressions for the low temperature residual 
entropy. We compute the static structure factor and find evidence for long range partial order for 
antiferromagnetic J2 , and short range magnetic order otherwise. The magnetic susceptibility crosses 
over from following a Curie- Weiss law at high temperatures to a low temperature Curie law whose 
slope clearly distinguishes ferromagnetic J2 from the J2 = case. We briefly comment on a recent 
report of the creation of a |th site depleted triangular lattice cobalt hydroxide oxalate. 

PACS numbers: 75.10.Hk, 75.40.Cx, 75.40.Mg 



I. INTRODUCTION 

The study of geometrically frustrated magnetic ma- 
terials provides a fruitful path to the discovery of ex- 
otic physics. In such materials, competing interactions 
between magnetic moments due to their localization on 
a highly symmetric lattice structure lead to large num- 
bers of quasi-degenerate states. Interesting physics re- 
sulting from such near degeneracies on frustrated lat- 
tices is thought to include: heavy fermion behavior in 
the spinel LiV^Oj^, evidence for the existence of mag- 
netic monopoles in the spin ice pyrochlores Dy2Ti20f^ 
and Ho2Ti20tH' ) non-Fermi liquid behaviour in the dis- 
torted kagome based FeCrAs^, spin liquid behavior in 
hyperkagome Na4lr 3 0j^ and triangular lattice based k- 
(BEDT-TTF) 2 Cu 2 (CN)J7 I , and unusual yet to be under- 
stood Kondo physics in the pyrochlore P^h^O^. 

While in traditional magnetic materials a two- 
dimensional net might arise in a layered lattice struc- 
ture, a growing body of work has focused on designing 
model systems to test models of strongly correlated elec- 
tron systems. The study of two-dimensional frustrated 
systems has be en au gmented by the creation of artificial 
spin ice systems^^, where people have shown that it is 
possible to lithographically etch ferromagnetic islands in 
patterns that can be easily modeled. From a theoretical 
perspective, the dipolar interactions crucial to such spin 
ice physics converge quickly in two dimensions, so that it 
is feasible to sum the dipolar interactions in real space. 
Recent work on triangular optical lattices^, magnetic 
colloids 12 and macroscopic spin ice^have generated fur- 
ther interest in two dimensional frustrated spin models. 

The Ising model on the two dimensional triangular net 
has been studied since the 1950'a^ and been shown to 
possess a large spin degeneracy as evidenced by a resid- 
ual ground state entropy per spin of ~ = 0.32cP^. The 



simplest frustrated depletion of this edge-shared bipar- 
tite net removes every fourth spin in a regular fashion, 
substituting a non-magnetic atom for every second atom 
along each of the basis vectors, to create the kagome net. 
The kagome net possesses a larger ground state degen- 
eracy for the Ising model than the triangular net due 
to the corner-sharing nature of its bonds, as shown by 
an enlarged residual entropy per spin of = 0.50183^. 
After many years of experimental searches, the kagome 
net has been experimentally realized both in quantum 
spin systemiP^lEi] an d artificial spin ice systems^!. The 
triangular kagome net is a j^th depleted triangular net, 
removing every fourth spin along the basis vectors and 
their neighboring hexagons. Once again, the antiferro- 
magnetic Ising model features a large residual entropy 
per spin, f = § ln(72) « 0.475185P The maple leaf 
nelP^is a =th depleted triangular net, removing a trian- 
gular lattice formed by advancing along one basis vec- 
tor by two units and a second by one unit in a manner 
similar to a knight in chess. The Ising model on this 
edge-sharing net, which has coordination number 5 in 
contrast to the triangular (6) and triangular kagome and 
kagome (4) coordinations, appears to have not yet been 
studied. Experimental candidates for both the triangular 
kagome — and maple leaPSHZS ne tg have been found. 

In this paper we treat the Ising model on the sorrel net, 
a l/9th site depleted and further 1 /7th bond depleted tri- 
angular corner-shared triangle net as described below in 
Sect. |TT] By means of Monte Carlo simulations, we find 
that the classical antiferromagnetic Ising model on this 
net has a finite residual entropy per spin in the thermo- 
dynamic limit of ■§ = 0.48185±0.00008, slightly greater 
than the Pauling entropy estimate for this net, indicating 
that the Ising model on this net is more highly frustrated 
than on the triangular and triangular kagome nets, even 
approaching that on the kagome net. We calculate the 
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FIG. 1: (Color online)(a) Selective removal of l/9th of the 
sites of the triangular net gives an new edge-shared net. The 
further removal of l/7th of the bonds (indicated by dashed 
lines labeled J2) results in the sorrel net shown in (b). There 
are 8 sites in the unit cell, six of which are 4-coordinated 
{a, f3, e, f , n, 1}, and 2 of which are six coordinated {7, 5}. (c) 
A photograph of oxalis regnelli triangularis, commonly known 
as wood sorrel with lines drawn atop it to guide the eye. (d) 
a common pattern on a doumbek. 



heat capacity, acceptance rate, magnetic susceptibility 
and static structure factor as a function of temperature 
to aid the experimental identification of candidate sys- 
tems. We investigate the stability of this model to the 
presence of magnetic interactions on the removed bonds 
of ferromagnetic or antiferromagnetic sign, and show that 



for small bond strengths the residual entropy in the spin 
system approaches that of the sorrel net to quite low 
temperatures prior to undergoing a secondary transition 
to a state still possessing a finite residual entropy in the 
T -> limit. 



II. A NEW TWO-DIMENSIONAL LATTICE 

In this paper we introduce a l/9th depleted triangu- 
lar net as shown in FigQa) of symmetry P6mm. One 
sees that the substitution of non-magnetic atoms (which 
themselves form a triangular lattice^) for every third 
atom along each of the basis vectors of the triangular 
lattice leads to an edge-shared net with sites 6-fold and 
5-fold coordinated. Amazingly, the symmetric selective 
removal of 1 /7 of the bonds on this lattice leads to a new 
highly frustrated corner-shared equilateral triangle two- 
dimensional net as shown in Fig. [ljb). This creates a 
net reminiscent of oxalis regnelli triangularis (commonly 
known as wood sorrel) leaves as pictured in FigjTJc) , and 
for this reason we will henceforth refer to this kagome- 
inspired net as the sorrel net. Fig. [T] (d) shows a pic- 
ture of a doumbek. The corner-shared triangle struc- 
ture pictured in Fig. [I] (b) is often found on the side of 
these drums, which are traditionally covered by mother- 
of-pearl, using simple Islamic tiling patterns. 

On the sorrel net the magnetic sites are four and six 
coordinated, such that one has two types of sites: one 
connecting two corner-shared triangles, another connect- 
ing three corner-shared triangles. To our knowledge this 
net has not yet been found in nature. It is our hope 
that our identification of this structure as an interesting 
frustrated lattice will encourage experimental groups to 
investigate materials naturally forming in this pattern, as 
has been the case with recent experimental searches for 
the maple leaf and kagome lattices. Indeed, one recent 
studjffl has managed to achieve a |th depleted triangu- 
lar lattice in a magnetic system, although, as we discuss 
below in Sec. |VI[ the antiferromagnetic order seen may 
have resulted from additional hexagonal layers in this 
complex material. 



III. MODEL 



We consider an antiferromagnetic Ising model on the 
sorrel lattice with the Hamiltonian, 



H= Ji^oto* + J 2 Yoio%, 

(ij) (ik) 



(1) 



where J\ > is an antiferromagnetic interaction between 
nearest neighbors on the corner-shared sorrel lattice, and 
Ji is an interaction between nearest neighbor sites that 
when coupled create an edge-shared net. In our simula- 
tions we allow Ji to vary in strength from —J\ to J\. 
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TABLE I: Sites, positions and nearest neighbors for atoms in 
the unit cell of the sorrel net. When not in the first unit cell 
neighbors are connected by basis vectors Ai = (a, 0, 0) and 
A2 = (§, jpjO). Ji > is an antiferromagnetic bond be- 
tween corner shared neighbors while J2 connects the remain- 
ing nearest neighbors and is treated in three limits:{ — Ji < 
J 2 < 0, Ji = 0,0 < J 2 < Ji}. The limit J 2 = J\ corresponds 
to the antiferromagnetic Ising model on the | th depleted tri- 
angular net. 



IV. METHOD 

We use the metropolis Monte Carlo algorithm to sim- 
ulate 8 x L x L nets subject to periodic boundary con- 
ditions. Lattice sizes have been varied from L = 1 to 
18 on the sorrel lattice with J2 =0 to ensure that the 
thermodynamic limit has been reached. At each temper- 
ature, prior to averaging, 2 x 10 5 to 10 6 Monte Carlo 
steps have been used to equlibrate the spin system, with 
each average then taken over a further 2001 to 20001 
Monte Carlo steps, where each step corresponds to on 
average attempting an update on each of the spin sites 
once. To simulate the annealing^ of an experimental 
material from high temperatures we have started with 
the temperature equal to 6000 J\, reducing the temper- 
ature by multiplying by a factor of 0.99 for successive 
measurements . 



A. Acceptance rate 

1. Single spin flip 

For the pure sorrel net (J2 = 0) we have carried out 
single spin flip Monte Carlo simulations and find a non- 
negligible acceptance rate for all L (rj 41.5% for L > 2) 
even in the zero temperature limit as shown in Fig. [2] 
(a). For small system sizes (L = 1, 2), we can find exact 
results for the zero temperature acceptance rates (0.425 
and 0.4148376 respectively). These are shown in the inset 
to Fig. [2] (a) to agree well with our simulations. The 
L = 2 acceptance rate appears to lie slightly lower than 
the acceptance rate for higher values of L, which appear 
to reach a constant value well before L — 15. 

When J 2 > the single spin flip process continues to 
admit flips on the 6 coordinated sites to the lowest tem- 
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FIG. 2: (Color online) (a) Acceptance rate versus tempera- 
ture for Ji — 1, J2 = 0. Representative lines for L — 1,2,9 
and L = 15 are shown. In the inset, we see that the zero 
temperature L=l acceptance rate is considerably larger than 
the acceptance rate for L — 2 and above. Exact ground state 
results for L = 1, 2 calculated as discussed in the text agree 
with the T — > limit of our simulations as shown by the 
horizontal lines in the inset and X's in the main figure, (b) 
Acceptance rate for J2 7^ [L = 4) with hybrid spin flip pro- 
cesses allowable as discussed in Sect. |IV A2| For J 2 > 0, as 
T — >0, twenty-five percent of all attempted single spin flip pro- 
cesses are accepted, corresponding to a magnetically ordered 
ground state with completely free spins at the 6-coordinated 
sites. For J 2 < single spin flip processes freeze out at low 
temperatures, while triple spin flip and double spin flip pro- 
cesses continue to allow the system access to all ground states 
at low temperatures. 



peratures, approaching a 25% acceptance rate. As the 
degrees of freedom associated to the spins at the remain- 
ing sites freeze out (as one expects from the detailed na- 
ture of the ground state configurations shown below, see 
Fig. [4]), the acceptance rate is not a monotonic function 
of temperature, dipping below 25% prior to recovering as 
T -> 0. 

When J 2 < 0, single spin flip simulations are insuffi- 
cient to allow the system access to its degenerate ground 
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state, and the acceptance rate dies quickly as T — > 0, in- 
dicating a lack of ergodicity as has been commonly found 
in the dipolar spin ice systems, where loop algorithms 
have commonly been introduced.^ 



2. Hybrid spin flip 

To address the concern that single spin flip processes 
alone may not allow our spin system to properly equili- 
brate for J 2 7^ 0, given that in the ground state the spins 
joined by the J 2 bonds want to align(anti-align) when 
J2 < (>)0, we have developed a hybrid spin flip algo- 
rithm. The addition of these flipping processes is in the 
spirit of loop algorithms 31 which have been deemed nec- 
essary in other frustrated magnetic systems to restore 
ergodicity to the system at low temperatures. To im- 
plement this hybrid spin flip algorithm, at each update 
attempt we randomly choose a site in the lattice. If this 
site belongs to a J 2 bond, we randomly choose to either 
flip the spin (single spin flip (ssf)), or to flip the pair of 
J2-bonded spins (double spin flip (dsf)). If the chosen 
site belongs to a 6 coordinated site, then we randomly 
choose to either flip the spin (ssf) or to flip the spin and 
one of its randomly choosen nearest J 2 bonded pairs of 
neighboring spins (triple spin flip (tsf)). Since 3/4 of 
all sites in the lattice belong to a J 2 pair, dsf processes 
are attempted on average during 3/8, tsf processes 1/8, 
and ssf processes 1/2 of the update attempts. With this 
hybrid code, one sees (in Fig. [2] (b)) that the total ac- 
ceptance rate remains finite to zero temperature for all 

J 2 ¥= 0. 

a. J 2 > As was found using ssf for an antifcr- 
romagnetic bond between the J 2 bonded sites, at the 
lowest temperatures only 1/4 of the spins remain active, 
and 1/4 of the ssf attempts within the hybrid spin flip 
are accepted. However, the presence of dsf and tsf pro- 
cesses does allow the temperature evolution of the ssf 
acceptance rate to become smooth and monotonic. This 
indicates that while the dsf and tsf processes effectively 
turn off once the system has achieved a partially or- 
dered ground state (as described below), their presence 
at higher temperatures allows the system to remain in 
equilibrium as this order occurs. 

b. J 2 < At low temperatures with a ferromagnetic 
bond between the J 2 bonded sites, ssf processes are no 
longer accepted. However tsf and dsf processes allow the 
system to access its degenerate ground states and move 
from one to another. 

c. J 2 — To test the validity of this approach, we 
have run hybrid spin flip processes for J 2 = 0. These give 
results which agree with those of the ssf code. 
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FIG. 3: (Color online) (a) Heat capacity vs. temperature for 
J'2 =0. For L > 2 heat capacities are very similar. For J2 7^ 
((b) and (c)) a sharp peak is seen below the J2 = peak. For 
J'2 > 0, (b), this peak sharpens with increasing system size. 
For J2 < 0, (c), both peaks show remarkably little variation 
with L. 



V. RESULTS 
A. Heat capacity 

As shown in Fig. [3] (a) , the heat capacity, 

r m < E2 >-< E > 2 
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FIG. 4: A phase diagram illustrating the robustness of the 
sorrel lattice entropy to variations in J2. For J2 > at low 
temperatures one sees completely unconstrained spins on the 
6 coordinated sites as illustrated. For J2 < below the transi- 
tion, satisfying the ferromagnetic J2 bonds leaves the system 
with considerable remaining degeneracy. 



for the pure sorrel net shows a broad peak indicat- 
ing the entrance of the system to a degenerate ground 
state which features two spins up(down) and one spin 
down(up) on each of the corner-shared triangles. Here 
E is the energy of a particular spin configuration on the 
lattice and T is the temperature. The heat capacity ex- 
hibits a mild dependence on system size, with the L = 1 
peak being noticeably larger than subsequent L. There 
is little variation in the heat capacity for L = 2 up to 
L = 18. In Fig. [3] (a) we show heat capacity versus tem- 
perature curves for L = 2 and L = 12 to illustrate that 
finite size effects are likely small. 
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FIG. 5: (Color online) Heat capacity per spin at J2 7^ minus 
the J2 = heat capacity per spin versus scaled temperature 
for L — 2. (a) For J2 > the peak seems to occur at the 
same scaled temperature, indicating a linear dependence of 
the transition as shown in Fig. [4] (b) For J2 < the crossover 
occurs at a lower temperature, and for small J2 maintains a 
linear J 2 dependence at the critical temperature. 



B. Phase diagram 



Turning on the coupling J2, we see a sharp peak in 
the heat capacity well below this broad feature for both 
signs of J2. As J 2 increases, this sharp feature moves 
to higher temperatures, broadens, and becomes less dis- 
tinct from the peak representing the crossover to a two 
up (down) and one down (up) spin state. In both cases, 
this secondary peak eliminates much of the entropy as- 
sociated to the pure sorrel net at temperatures which 
increase with the magnitude of J2- As seen in Fig. [3] 

(b) , when J2 > 0, as L increases the second peak in the 
heat capacity sharpens as one might expect from a phase 
transition to a long range ordered state. In strong con- 
trast to this behavior, when J2 < as seen in Fig. [3] 

(c) , the second peak remains essentially unchanged as L 
increases, indicating a crossover to a short-range ordered 
state. 



In Fig. [4] we sketch the phase diagram for the Ising 
model on the generalized sorrel net. In this plot, the 
dashed line represents the crossover from local disorder 
to each Ji triangle having a net energy —J\. That is, be- 
low this line each triangle of the pure sorrel lattice adopts 
a spin structure with two spins up (down) and one spin 
down (up). From the heat capacity, the temperature as- 
sociated to this crossover corresponds to the maximum 
seen in the broad first peak of C v versus T, which ap- 
pears not to change much as J2 is turned on. As the 
temperature is further decreased, for J 2 > the spin 
system undergoes a transition to a long range partially 
ordered state, as pictured in Fig. [4] In this state the 
6-coordinated spin sites are completely free Ising spins, 
while all other nearest neighbor spin pairs feature an- 
tiparallel spins. 
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For J2 < 0, as the temperature decreases there is 
a (size independent) crossover to a short-range ordered 
state, where across each J 2 bond the spins are parallel. 
Every pair of parallel spins along J2 bonds is possible, 
with the 6-coordinated spin site adopting a unique spin 
configuration antiparallel to the dominant spin configu- 
ration of its 3 nearest neighbor pairs. One sees that such 
states are connected by the above-described dsf and tsf 
updates, and not by ssf updates. 

At J2 = 0, the spin system remains disordered to zero 
temperature classically. As J2 turns on, the peak cor- 
responding to the transition/crossover to a long range 
ordered/short range ordered state occurs at proportion- 
ally larger temperatures, as indicated by the straight lines 
in our sketched phase diagram. If the high temperature 
crossover and its associated temperature dependent heat 
capacity is truly independent of J2, then by subtracting 
out the J2 = peak one would be left with the second 
peak alone. Further, if one then scaled the temperature 
of several of these curves by their value of J 2 , their peaks 
would be expected to coincide in the case that the transi- 
tion temperatures scale linearly with J2. Such plots (for 
L = 2) are shown in Fig. [5] (a) (for J 2 > 0) and (b) (for 
J 2 < 0). One sees that the temperature corresponding to 
the peaks of these curves does show a remarkable coinci- 
dence in both cases. In fact, only the peak of J2 = —J\ 
appears to lie marginally below this linear temperature 
dependence. 



C. Entropy 

At very high temperature one expects there to be es- 
sentially no correlations between the Ising spins, or 2 N 
equally weighted spin configurations for an entropy per 
spin of ln(2). As one lowers the temperature, the resid- 
ual entropy in the system decreases. In ordered materials 
below the ordering transition temperature the entropy 
quickly approaches zero. In magnetically frustrated spin 
systems entropy can remain to very low temperatures. 
As such it is interesting to investigate the behavior of 
the residual entropy in the system, 



S(T) = S(T = 00) 



' Cy(T) 

T 



dT, 



(3) 



and the zero temperature limit S(T = 0) of this quantity. 



1. Pauling estimate 

When J2 = 0, if one makes the assumption that spins 
beyond their nearest neighbor are uncorrelated, then 
each spin site has two degrees of freedom constrained by 
the fraction of the possible spin configurations on each 
triangle which belong to the ground state, which for a 
triangle-based lattice is |. On the sorrel net each spin 
belongs to either two or three triangles, so that there 
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FIG. 6: (Color online) Residual entropy (as T — > 0) versus 
system size (L) is roughly constant beyond L — 2. Shown 
are data from two different equilibration times with the same 
initial seeds and the effect on the heat capacity of a single 
spurious low temperature point in the heat capacity. Simu- 
lations agree well with exact results at L = 2 and L = 1, 
which is shown on a different scale. (Inset) Residual entropy 
versus for L > 4 is used to extract a residual entropy 



= 0.48185 ± 0.00008 in the thermodynamic limit. 



are N t = ^jf triangles in total. The Pauling degener- 
acy of the ground state is then 2^(1)^* = 2 N Q)™, 

hence the Pauling estimate for the entropy per spin is 

S ■ ■ - 

N 



= ln(2(f)l) = iln(f)« 0.477386 



2. Exact 

When the number of spins is small, it is numerically 
possible to count the number of ground states that satisfy 
the ground state conditions for J 2 = 0. For any spin 
configuration in the ground state two of the spins on any 
given triangle must be up (down) while the remaining spin 
is down(up). So doing, for L=l we find the ground state 
to have 40 distinct members while for L=2 it has 4921350 
distinct states. This indicates that the exact result for 



the residual entropy of the pure sorrel net is 



JV 



0.46111 



for L=l and § = 0.481534 for L=2, both numbers which 
fall within our error bars for the residual entropy as seen 
in Fig. [6] 

Adding J 2 ^ 0, the available ground states addition- 
ally must satisfy parallel(antiparallel) spins across the J2 
bond when J2 < (>)0. This means that each ground 
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FIG. 8: (Color online) The entropy remaining in the system 
per spin as a function of temperature for L = 4. At high 
temperatures, all curves asymptotically approach ln(2). For 
J2 = 0, the residual entropy smoothly decreases until reaching 
a constant low temperature value -jl = 0.48173±0. 00005. For 
small J2, the entropy remaining follows the J2 = curve to 
quite low temperatures prior to sharply downturning to the 
exact T — residual entropies shown in the text. 



have exact results for the residual entropy. 
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FIG. 7: (Color online) As discussed in the text, we have found 
exact results for the ground state residual entropy for J2 7^ 0. 
These agree well with MC simulations for (a) J2 > and (b) 
J 2 < 0. 



state is a subset of the sorrel lattice states. For L = 1 
this reduces the number of degenerate ground states to 
just 8 independent of the sign of J 2 . For L — 2 when 
J 2 > only 512 of the ground states satisfy the con- 
dition, while for J 2 < there are 4096 ground states. 
We see that these numerically extracted exact values 
are the L = 1 and L = 2 limits of general exact for- 
mulae for the J 2 7^ entropy. For J 2 > all the 
ground states feature a completely free spin on the 6 
coordinated sites, and only one of two global spin con- 
figurations for the remaining spin sites, indicating that 
we have 2 2L +1 degenerate ground states, and that the 
entropy per spin should therefore weakly depend on L: 
s = ^l n (2^ 2 +i) = (i + g^)ln(2). For J 2 < 0, the 
ground states feature parallel spins on each of the J 2 
bonds, meaning that each 6-coordinated site has a defi- 
nite orientation antiparallel to the dominant spin type of 
its three neighboring J 2 bonded pairs of spins. However, 
every combination of J 2 bond parallel spins is possible, 
so there are 2 3i ground states, for a ground state en- 
tropy per spin of = 31 'g 2 - ) . Hence for all J 2 7^ we 



3. Monte Carlo 

As shown in Fig. [6j the residual entropy per spin calcu- 
lated by Monte Carlo simulations'-? agrees well with the 
exact results at L — 1 and L — 2 for the pure sorrel net. 
Fig. [6] additionally presents data using both the ssf and 
hsf codes. We see that simulations from both codes agree 
with the exact results and give equally valid results at 
larger L. We further note that for L = 2 and beyond, the 
residual entropy per spin for this case (J 2 = 0) appears 
to have reached a constant value (0.48185±0.00008p21 as 
a function of system size, indicating that the thermody- 
namic limit has been reached well before our maximum 
system size of L — 18, which corresponds to 2592 spins. 

As expected from our exact results, much of this resid- 
ual entropy is quenched at very low temperatures when 
J 2 7^ 0. Our simulations indicate that the residual en- 
tropy per spin for J 2 < remains constant as a func- 
tion of L. All calculated residual entropies for J 2 < 
agree within error of the exact results, except for the 
case J 2 = — J 1 and L = 4, where a marginally lower en- 
tropy is found. The reason for this small discrepancy has 
not been discovered as yet. To calculate the errors for 
these results we numerically integrated four (or more at 
small L) heat capacities divided by temperature to find 
entropies versus temperature. We then found the average 
entropy at each temperature and its standard deviation. 
As such, our stated errors may be a little small as they 
do not explicitly account for errors in the numerical in- 
tegration procedure. 



T = 5.07 J 




FIG. 9: For L = 6, an evaluation of the static structure 
factor at the indicated temperatures in reciprocal space for 
J'2 = 0, contour increments of 0.1 for (a) and 0.2 for (b). 
The squished hexagon at the center represents the first Bril- 
louin zone, which has been repeated to tile the reciprocal 
space. The center of each hexagon indicates the number of 
reciprocal lattice vector steps to the centre from each side, 
the top center being (0,-^^). One sees diffuse weight about 
q = ( ^ , 0) and the 5 equivalent Ce rotations about the origin 
at all temperatures. 
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FIG. 10: (Color online) The (L = 6) static structure factor 
for J 2 = at T = 0.251Ji and 10.05 Ji along selected high 
symmetry directions. Note that at both temperatures the 
weight remains symmetric about the origin to about q = . 
Along (0, q) the structure factor then repeats (note that this 
distance is 1.5 reciprocal lattice vectors). Diffuse maxima 
occur about the points q = (±— ,0) which are considerably 
less distinct at higher temperatures. 



D. Spin-spin correlations 



For J 2 > the residual entropy per spin is consistent 
with the exact results with all cases lying within two 
standard deviations.^ 

In Fig. [8] we plot the entropy versus temperature for 
various values of J 2 at L — 4. We see that at high tem- 
peratures as J2 is turned on from zero, all the curves 
follow the J2 = entropy. As the temperature continues 
to decrease, the larger magnitude J2 simulations begin 
to substantially deviate from the J2 = entropy in a 
manner consistent with the phase diagram as extracted 
from the heat capacity in Fig. |lj When | J%\ is small, this 
deviation does not occur until quite low temperatures, 
and when it does occur, the transition is quite sharp to a 
state with the ground state entropy. If a system can be 
made for which the Ising model on the sorrel lattice with 
J 2 « J\ is a good approximation, one would expect to 
see a plateau in the residual entropy of the pure sorrel 
lattice over a certain temperature range, as shown in the 
inset to Fig. [8j 



The static structure factor provides a measure of the 
strength of the spin-spin correlations and their momen- 
tum dependence. Where large single crystals of a ma- 
terial become available, the structure factor (convolved 
with an atomic form factor) allows a direct comparison 
between neutron scattering measurements and theoreti- 
cal predictions. Materials which order magnetically show 
the development of sharp, resolution-limited magnetic 
Bragg peaks at low temperatures. In contrast, many frus- 
trated magnetic systems show evidence of short range or- 
der in the form of disperse magnetic peaks concentrated 
around certain areas of reciprocal space. Powder neutron 
scattering measurements average the angular dependence 
of such graphs, showing weight at finite wavenumbers. 

Within Monte Carlo, the static structure factor is de- 
fined as, 

S ^ T ) ^ ^E <S ^( T ) ' «J(T)> e**^-^, (4) 
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T =5.07 J ! 




FIG. 11: For L = 6 an evaluation of the static structure 
factor at the indicated temperatures for J2 = Ji in reciprocal 
space. Contour increments of 0.2 in (a) precede the resolution- 
limited Bragg peaks of (b). The center of each hexagon in- 
dicates the number of reciprocal lattice vector steps to any 
point. It is four reciprocal lattice vectors to the center from 
each side, the top center being (0, ^^)- 



where N is the number of sites in the lattice N = 8L 2 , 
the summation runs over all pairs of spins on the lattice, 
the braces indicate the average over several MCS (20001 
here) of the dot product between spins at sites i and j, 
and fl — fj is the vector between spin locations at sites 
i and j. For the sorrel lattice we have a hexagonal lat- 
tice with an eight site basis and reciprocal lattice vectors 

bi = ?f (l> -^) and b 2 = ^(0, 1), where a is the di- 
rect space lattice constant. For a finite 8 x L x L lattice, 
we can evaluate the structure factor at q = jb\ + jbi, 
where / and J are integers. Taking / and J to run inde- 
pendently from {— 4L,4L} allows us to gain a picture of 
the magnetic correlations well beyond the first Brillouin 
zone (shown as the central squished hexagon in Fig. KM. 
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FIG. 12: For J2 = Ji the temperature dependence of the 
peak of the structure factor at q = (— , 0) and equivalent 
points. By T — 0.251Ji it has almost reached the ideal or- 
dered value as discussed in the text. 



along the (g,0) direction at q = ±— and the equivalent 
q = ±— and symmetry equivalent (Cq rotation) points. 
This corresponds to a linear combination of reciprocal 
lattice vectors, for instance, 2&i +62- As the temperature 
is lowered, the relative intensity of this peak increases 
but the weight remains diffusely distributed around this 
point. 



2. J 2 > 

For J2 at high temperatures the structure fac- 
tor appears quite similar to the J2 = case, as seen 
in Fig. 11 (a) and Fig. 



13 (a). However, for J 2 > 0, 



as the temperature is lowered one sees the rapid shift- 
ing of weight to the maximal direction-which becomes 
an ordering wavevector as shown in Fig. 11 (b). In Fig. 
|12| we present the temperature dependence of the struc- 
ture factor at this wavevector. One sees that it saturates 
at a value of approximately 121.70 by the temperature 
T = 0.25Ji. In this partially ordered state we have ar- 
gued that the spins at the 6-coordinated sites are com- 
pletely free to fluctuate while the remaining spins form a 
magnetically ordered state. If one simply assumes a spin 
structure with no moments on the 6-coordinate sites and 
10 A disperse maximum is seen antiferromagnetically ordered classical spins on the re- 



1. J 2 = 

When J2 = 0, the structure factor weight is lowest in 
the vicinity of q = 0, and isotropic at fixed \q\ as seen 
in Fig. [9] and Fig 
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maining sites, one finds a structure factor of 121.50. It 
is interesting to note that this is not the full moment 
| • 288 = 216 corresponding to an alignment along a par- 
ticular axis of the spin structure, as one is used to from 
both ordered ferromagnets and antifcrromagncts. In the 
partially ordered antifcrromagnetic state although it is a 
two-sublattice structure, one sublattice gains a weight- 
ing of \ on its moments due to the geometry of the lat- 
tice in its structure factor, lowering the total moment. 
Adding back the self correlation of free (unit magnitude) 
spins gains an additional 0.25 contribution to this ideal- 
ized structure factor, suggesting that by T = 0.25Ji the 
ordering state may have nearly reached the ideal partially 
ordered structure. 




3. J 2 < 

When J2 < 0, the low temperature spin correlations 
(Fig. 13 (b)) show once more (like when J2 = 0) a broad 
feature about q = (±47r, 0), q = (±87r, 0) and 6- fold sym- 
metry related points. This feature sharpens as the tem- 
perature lowers, as can be seen in Fig. [l4j Additionally, 
the contours appear less circular, more hexagonal than 
the J-2 = case. 




E. Magnetic susceptibility 

The static magnetic susceptibility is one of the first 
measurements made on new magnetic materials. It is 
calculated as, 



X(T) 



S(q = 6,T) 1 



T 



NT 



2<«J(T).sJ(T)>. (5) 



1,3 



Plots of the inverse magnetic susceptibility versus tem- 
perature are shown in Fig. |15| At high temperatures 
(Fig. 15 (a)), the inverse susceptibility shows a Curie- 
Weiss form with a large negative Curie- Weiss constant 
indicating the dominance of antiferromagnetic couplings 
between nearest neighboring sites. A characteristic fea- 
ture of frustrated magnets is that their magnetic suscep- 
tibility shows no signs of magnetic order to well below 
the Curie- Weiss temperature. Unlike ordering antiferro- 
magnets which exhibit non-monotonic susceptibilities at 
the Neel temperature, for —J% < J2 < Ji, the inverse 
susceptibility monotonically decreases with the tempera- 
ture. It is not uncommon for frustrated materials to have 
defects or uncorrelated "orphan spins" , such that the sus- 
ceptibility diverges more quickly at low temperatures, as 
a Curie law. In this system we see an intrinsic property 
of Ising spins on the sorrel lattice is a low temperature 
Curie tail. 



FIG. 13: For L = 6 an evaluation of the static structure factor 
at indicated temperatures for J2 = — Jj in reciprocal space, 
contour increments of 0.1 (a) and 0.2 (b). The center of each 
hexagon indicates the number of reciprocal lattice vector steps 
to any point. It is four reciprocal lattice vectors to the center 
from each side, the top center being (0, ^Jg^)- (a) At T = 
5.07Ji, the highest weight lies at (0, ±^|^) and equivalent 

positions, but considerable weig ht lies at q = (±v>°) which 
grows at low T to form diffuse peaks. 



1. J2 = 

When J 2 = 0, the inverse susceptibility crosses over 
from a high temperature Curie- Weiss law to a low tem- 
perature Curie law. For L = 2 we find x 1 = (3.9134 ± 
0.0007)T. This likely indicates that a substantial fraction 
of the spins at the 6-coordinated sites are freely flipping 
at low temperatures. 



2. J2 > 

For J2 > 0, this is certainly true. At the lowest tem- 
peratures (see Fig. [l5"|h)) all J 2 > curves that we have 
investigated approach the same Curie law. For J2 = 1.0, 
and L = 4 we find this to be x _1 = (4.000 ± 0.001)T. 
This is consistent with | of the spins being completely 
free Ising spins, while the remaining spins are ordered. 
Interestingly, one does see that at temperatures of order 
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FIG. 14: The structure factor at low and high temperatures 
is now roughly symmetric only from 



4=- to -^S- and is less 



rounded than at J2 = (Fig. 10 1. A diffuse maxima occurs 
at q = (± 0) at low temperatures but washes out at higher 
temperatures. 



—dew = (5.97 ± 0.02) Ji, the inverse susceptibility begins 
to deviate upwards from its Curie- Weiss fit, as one might 
expect from an ordering antiferromagnet. However, as it 
begins to order, it frees the 6-coordinated spins, causing 
X 1 to eventually downturn at low temperatures once 
more. As shown in Fig. 15 (a), there is surprisingly little 
system size-dependence to the magnetic susceptibility, so 
we have chosen to show L = 2 plots in Fig. 15 (b). 




J = -0.05 1 



u 0.01 0.02 0.03 0.04 
T[Jj] 



3. J 2 < 



The low temperature susceptibility for J2 < follows a 
distinctly different Curie law. For L — 2 it is not hard to 
write out all possible ground state spin structures consis- 
tent with every J 2 bond having aligned spins. So doing, 
one finds total spins ranging from S = to S = 16 
and their degeneracies as shown in Table [TTJ Writing the 
magnetization in the limiting case B —¥ in terms of 
the magnetic moment per spin state and the partition 
function allows us to evaluate the magnetic susceptibil- 
ity and find an exact result, \l=2 = fi T ~ 2.2857T. 
Our Monte Carlo results^ at L = 2 agree with this 
exact calculation for J2 = —J\ at low temperatures, 
X^HT) = (2.28689±0.00059)T-(7±30)xl0~ 5 . The low 
temperature asymptotes for J < for higher L appear 
to agree with this result as well. 



FIG. 15: (Color online) Inverse magnetic susceptibility ver- 
sus temperature, (a) A high temperature Curie- Weiss inter- 
cept gets increasingly negative as J2 increases. Shown are 
L — 2 and L = 4 for all curves, and additionally L = 6 for 
J2 = —J\. There is very little system size-dependence of any 
part of the susceptibility beyond L — 2. (b) Variation of J2 
for L — 2. Curves asymptote at low temperatures towards 
different Curie laws depending on the sign of J2, approaching 
less quickly the weaker J2 is. (c) The Curie law for J2 = 
slightly differs from that for J2 > 0. 



VI. DISCUSSION 
A. Non-magnetic depletion 

In this work we have introduced the sorrel net and 
provided the first study of frustrated magnetism on this 
depleted triangular lattice. The regular substitution of 
non-magnetic atoms into a frustrated magnetic sublattice 
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Total spin 


Degeneracy 


±16 


1 


±12 


12 


±10 


24 


±8 


96 


±6 


200 


±4 


516 


±2 


672 





1054 



TABLE II: Degeneracies of the ground state spin configura- 
tions satisfying parallel spins on every pair of J2 bonded spins 
as a function of the total spin of the L — 2 spin configuration. 



lattice ^jth depletion of a triangular lattice, where each 
point of the triangular lattice removes 7 spins. 

The natural extension of the doping of the triangular 
lattice from the kagome lattice (which removes every sec- 
ond spin along the reciprocal lattice vectors), to the sor- 
rel lattice (which removes every third spin), would lead 
to the removal of every fourth spin along the reciprocal 
lattice vectors for a ygth doping. Interestingly, for such 
a lattice it is not possible to remove bonds to create a 
new corner-sharing triangle lattice. However, the selec- 
tive removal of edge-sharing bonds analogous to setting 
J2 = on the sorrel lattice, would create an interesting 
potentially frustrated lattice-featuring a star of David 
decoration inside each empty hexagon of the triangular 
kagome lattice. 



of a material is a useful way of predicting and perhaps 
generating new candidate exotic materials. One such ma- 
terial is the quantum spin liquid candidate Na4lr30g^, 
which features s = ^ Ir atoms at three of the four cor- 
ners of the corner-shared tetrahedral lattice (common 
to the pyrochlore and B-site spinel geometries). The Ir 
atoms in this material are arranged to create a three- 
dimensional corner-shared triangle hyperkagome lattice. 
Non-magnetic Na atoms sit at the fourth corner of each 
tetrahedron, and whether because of substantial size 
differences^ between Na + and Ir 4+ ions or electrostatic 
interactions with surrounding atoms, the Na ions spread 
throughout the lattice in a regular pattern. 



B. Evidence for depleted triangular lattices 

In this context, it is interesting to note that tri- 
angular lattice depletions of the triangular lattice not 
only yield the (jjjth doped) sorrel lattice presented 
here, but also a number of frustrated lattices for 
which experimental candidates have already been found. 
An impressive collection of pictures of possible trigo- 
nally symmetric lattices derived by periodic site de- 
pletion is presented in the electronic supplementary 
information to Ref. [TJ Most notably, the much 
studied kagome lattice, a \ doped triangular lattice, 
has recently been realized experimentally in the ma- 
terials volborthite^, herbertsmithite^, kapellasite^, 
haydeeite^Sl, vesignieite^, and Cs 2 Cu 3 SnF 1 j22l. Many of 
these minerals are Cu based as experimentalists search 
for low spin quantum spin liquid candidates. Addition- 
ally, artificial spin ice on the kagome lattice was first 
produced by Qi et aZ.E2l 

The hth. depleted triangular lattice, the maple leaf lat- 
tice, was introduced in 1995 and has candidate materials: 
M rE [Fe(02CCH2)2NCH 2 P0 3 ] 6 • nH 2 where x = 11 if M 
is Na or K, and x = 10 if M is Rb 26 and spangolite! 211221 . 
Even the triangular kagome lattice realized by the Cu 
atoms^in CugX 2 (2-carboxypentonic acid)6-xrl 2 0, where 
X is F, CI or Br, studied by Loh et a/P2l is a triangular 



C. A |th depleted triangular lattice 

Is it possible to create a regular |th site depleted 
triangular magnetic lattice in an insulator? We were 
pleasantly surprised to learn that this has indeed re- 
cently been achieved in a new cobalt hydroxide oxalate 
Coi 2 (OH)i 8 (C2O4)3(C 4 N2Hi ) by Keene et oP. In this 
material the planar Co atoms are octahedrally oxygen 
coordinated, with | of the octahedra empty. Unfortu- 
nately, above and below the empty (in-plane) octahedra 
are Co atoms. Three oxygen atoms from the tetrahedra 
forming the base of a tetrahedral coordination of these 
Co atoms with nitrogen providing the apical atom of the 
tetrahedra. It is not immediately clear what effect cou- 
plings to such tetrahedral Co atoms might create, al- 
though even a weak coupling between the tetrahedral Co 
atoms might make a strong (ferromagnetic) perturbation 
on the hexagon of Co atoms nearest to them. Addition- 
ally, this layered material has a 2 dimensional honey- 
comb network of oxalate mediated couplings between Co 
ions. Magnetically, this material shows interesting, al- 
though perhaps not frustrated behavior, apparently an- 
tifcrromagnetically ordering at T c — 23.5 ± 0.5K despite 
an antiferromagnetic Weiss constant of only about 3.6 
K. Given that the author^ expect an antiferromagnetic 
coupling in the honeycomb layer of 18±3K from previous 
worlP3 it seems natural to ask whether or not the bonds 
of the sorrel layers in this material indeed have antiferro- 
magnetic couplings, or are contributing to lowering the 
effective Weiss constant. 

One expects that the predominant interaction between 
Co atoms in the plane is via superexchange through the 
joining oxygen sites, and that as the Co-O-Co bond varies 
from 180° to 90° at some angle (roughly 120° in the 
cuprates) ferromagnetic exchange will be favored over an- 
tiferromagnetic exchange. Of the four in plane bonds pre- 
sented by the authors, {J7i, J2, J3, Ja\ involving planar 
Co atoms, our Ji corresponds to the assumption of equal 
bond strengths on their J\ and J2, our J2 corresponds 
to their J% and their J± is a new bond to the tetrag- 
onal Co atoms. While it would be interesting to carry 
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out a detailed calculation of these exchange integrals in 
future work, naively from bond angle considerations^ 
one would say that J± is likely to be weakly antifer- 
romagnetic, J\ and J% may well be weakly ferromag- 
netic, and J-$ is likely to be weakly ferromagnetic. As 
such, one might expect that this system falls in a regime 
not covered in this work, of ferromagnetic J\ and J 2- 
It would certainly also be interesting if a local probe 
such as oxygen NMR could be used to determine the 
signs of the magnetic couplings within the plane. If, on 
the other hand, it were shown that this material pos- 
sesses an antiferromagnetic J\, one might hope to ex- 
ploit the differences between the atomic sizes of octahe- 
drally and tetrahedrally coordinated Co 2+ ions to design 
a non-magnetic substitution. For example, according to 
ShannoiPS, tetrahedrally coordinated Mg 2+ is very close 
in size to tetrahedrally coordinated Co 2+ , while the oc- 
tahedral coordinated ions of both species are of quite 
different sizes, suggesting that if it were possible to sub- 
stitute Mg atoms for Co atoms, the Mg atoms might 
preferentially enter the tetrahedral site, removing some 
complicating non-frustrated magnetic physics. 

D. Designing for antiferromagnetic superexchange 

While we view the solution of the Ising model on the 
sorrel lattice as a simple test whether or not antiferro- 
magnetism would be frustrated on this net, it is inter- 
esting to speculate on how this, or related models might 
be realized in an experimental system. In particular, one 
might wonder how magnetic exchange processes could 
lead to an antiferromagnetic J\ and a weak, non-existent, 
or ferromagnetic J2, even if the magnetic atoms of a 
crystal could be arranged into the ^th doped triangu- 
lar net. One could imagine^ superexchange pathways 
via say oxygen atoms living at the midpoints of the J\ 
bonds leading to a strong antiferromagnetic coupling on 
these bonds alone. Indeed even oxygen atoms lying at the 
center of the J\ triangles would likely lead to weakly an- 
tiferromagnetic Ji bonds. Such interactions might be ex- 
pected to produce at least roughly equal strength bonds 
around both the 6-coordinated and 4-coordinated spin 
sites, as we have assumed in this work. To addition- 
ally realize Ising spins would likely require strong crystal 
field anisotropies in the out of plane direction, leading to 
a magnetically easy axis. 

E. A two-dimensional structure 

That the sorrel net is two-dimensional opens a host 
of possible realizations for further study. For example, 
recent advances in the study of spinless bosons trapped 
by optical lattices have allowed the simulation of models 
of frustrated magnetism on the triangular lattice with 
couplings of various magnitudes!^ Could such studies 
be extended to the doped triangular lattices? From soft 



condensed matter, the study of magnetic colloids inter- 
acting via the dipolar interaction offers the possibility of 
real time manipulation of a lattice using optical tweez- 
ers to perhaps create a mesoscopic frustrated lattice of 
choice .121 

While it may be currently unrealistic^ to expect the 
creation of a direct artificial realization of a dipolar Ising 
model on the sorrel lattice with moments pointing out of 
the plane, the construction of a coplanar artificial spin 
ice structure should be feasible. As for (lithographically 
etched) realizations of artificial spin ice on the kagome 
lattice, a permalloy island could be to chosen to point 
locally along the axis of symmetry of the triangle. While 
such islands could join the centers of the triangles about 
the 4-coordinated spins sites, a new and interesting fea- 
ture for study would arise at the 6-coordinated spin sites. 
Three overlapping ferromagnetic islands (with a Cq sym- 
metry) would mimic a multiorbital local Ising spin at 
these vertices, producing a novel interplay between short 
range contact terms (Hund's coupling and ferromagnetic 
interactions) with the long range dipolar interaction. We 
have begurP21 the study of ferromagnetic interactions be- 
tween such local Ising spins in the presence of Hund's 
coupling terms and find, at least to this approximation, 
that degeneracy remains in the ground state of the sorrel 
lattice. Unfortunately, this multiorbital nature appears 
to prevent the creation of a macroscopic spin ice along 
the lines of that created by Mellado et aZpS 

The question of what quantum magnetism will do on 
the sorrel lattice is left for future work. When considering 
frustrated quantum spin systems, perhaps the simplest 
model is that of the quantum transverse Ising mode^, 
the study of which is considerably simplified for low di- 
mensional lattices. We^ have begun to investigate clas- 
sical dimer models on the sorrel lattice using the Pfaffian 
approach. The hope is that as on the triangular lattice, 
the classical and quantum dimer models may be related, 
and there may be a mapping^ to the quantum transverse 
Ising model. 



VII. CONCLUSIONS 

We have studied the antiferromagnetic Ising model 
on a new corner-shared equilateral triangle net, find- 
ing a large finite residual entropy j| =0.48185±0. 00008, 
slightly larger than the corresponding Pauling entropy 
jf = \ ln(") « 0.477386. We have shown that the spin- 
spin correlations remain disordered to the lowest temper- 
ature, with the structure factor showing broad disperse 
peaks about q = ( — ,0) and C§ rotations. The magnetic 
susceptibility indicates that strong antiferromagnetic cor- 
relations at high temperatures give way to short range or- 
dering at low temperatures, with a Curie tail with slope 
close to that one would expect were \ of the Ising spins 
completely free. While we have not proposed an exact so- 
lution, our results are remarkably independent of system 
size, indicating that the thermodynamic limit is quickly 
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reached. 

To this corner-shared triangle sorrel net we have 
added edge-sharing bonds (J 2 ) which are likely to remain 
present in physical realizations of the doped triangular 
lattice, and investigated the phase diagram produced by 
varying the magnitude and sign of this coupling constant. 

For antiferromagnetic edge-sharing couplings (J 2 > 0), 
we have shown that at low temperatures our spin system 
adopts a partially ordered state with Ising spins at the 
6-coordinated sites completely free, and all other neigh- 
boring spins antifcrromagnetically correlated, resulting 
in sharp magnetic Bragg peaks at q — ( , 0) and Cq 
rotations. This allows us to write exact solutions for the 
ground state entropy % = {\ + 8Z2) m (2) and the low 
temperature Curie tail of the susceptibility, results which 
are supported numerically. 

For ferromagnetic edge-sharing couplings (J 2 < 0), we 
have shown that at sufficiently low temperatures all com- 
binations of ferromagnetically aligned spins across edge- 
sharing (J 2 ) bonds are equally realized, resulting in dif- 



fuse magnetic scattering about q = ( — ,0) and Cq rota- 
tions, similar to the pure sorrel net case. In contrast to 
this pure case, we have shown that the low temperature 
susceptibility shows a considerably stronger Curie tail, 
which we have found to be consistent with exact L = 2 
results. Knowing the ground state, we have found an ex- 
act expression for the residual entropy ~ 3 1 of the 
spin system which is supported numerically. 
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